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ABSTRACT 



We present an atlas of Spitzer Space Telescope Infrared Spectrograph (IRS) 
spectra of highly luminous, compact mid-infrared sources in the Large Magellanic 
Cloud. Sources were selected on the basis of infrared colors and 8 /xm (MSX) 
fluxes indicative of highly evolved, intermediate- to high-mass stars with current 
or recent mass loss at large rates. We determine the chemistry of the circumstellar 
envelope from the mid-IR continuum and spectral features and classify the spec- 
tral types of the stars. In the sample of 60 sources, we find 21 Red Supergiants 
(RSGs), 16 C-rich Asymptotic Giant Branch (AGB) stars, 11 H II regions, 4 likely 
0-rich AGB stars, 4 Galactic 0-rich AGB stars, 2 OH/IR stars, and 2 B[e] super- 
giants with peculiar IR spectra. We find that the overwhelming majority of the 
sample AGB stars (with typical IR luminosities ~ 10^ Lq) have C-rich envelopes, 
while the 0-rich objects are predominantly luminous RSGs with Lir~ 10^ Lq. 
For both classes of evolved star (C-rich AGB stars and RSGs), we use the near- 
to mid-infrared spectral energy distributions to determine mean bolometric cor- 
rections to the stellar K-band flux densities. For carbon stars, the bolometric 
corrections depend on the infrared color, whereas for RSGs, the bolometric cor- 
rection is independent of IR color. Our results reveal that objects previously 
classified as PNe on the basis of IR colors are in fact compact H II regions with 
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very red IRS spectra that include strong atomic recombination lines and poly- 
cyclic aromatic hydrocarbon (PAH) emission features. We demonstrate that the 
IRS spectral classes in our sample separate clearly in infrared color-color dia- 
grams that use combinations of 2MASS data and synthetic IRAC / MIPS fluxes 
derived from the IRS spectra. On this basis, we suggest diagnostics to identify 
and classify, with high confidence levels, IR-luminous evolved stars and compact 
H II regions in nearby galaxies using Spitzer and near-infrared photometry. 

Subject headings: atlases — stars: AGB and post-AGB — Magellanic Clouds — 
infrared: stars — stars: mass loss — circumstellar matter 



1. INTRODUCTION 

While studies in visible light, particularly spectroscopy, have been key to our under- 
standing the nature and evolution of stars during most of their life cycles, the early and 
late stages of stellar evolution are best studied in the infrared. This is because, during these 
stages, the stars are embedded in regions of dust and gas. The late stages of stars of interme- 
diate initial mass (1-8 M©) are characterized by very high luminosities and obscuration by 
dusty, expanding circumstellar envelopes (see review by van Winckel 2003). Circumstellar 
dust envelopes absorb photospheric optical and ultraviolet emission and re-radiate it in the 
mid- to far- infrared (MFIR). Thus evolved intermediate-mass stars are luminous in the in- 
frared (IR) but may be highly obscured at optical wavelengths. Current theories of the later 
stages of the evolution of mass-losing stars therefore have been informed by IRAS, ISO, and 
2MASS studies of intermediate-mass evolved stars in the Milky Way Galaxy and Magellanic 
Clouds. 

Understanding the envelope chemistry of optically-obscured post-main sequence (MS) 
stars through IR studies is important for our overall view of galaxy evolution, since these 
stars dominate the rate of return of stellar processed material to the ISM (Jura & Kleinmann 
1989; Barnbaum et al. 1991; Le Bertre et al. 2001). Optical surveys (e.g., Groenewegen 1999 
and refs. therein) can miss these objects because they are too faint to be detected. Dusty 
circumstellar shells in evolved stars are produced as the star loses mass via a combination 
of pulsations and radiation pressure on dust grains newly formed via the pulsation process. 
Once accelerated, these grains strip the (predominantly molecular) gas away from the stellar 
envelope. While on the Asymptotic Giant Branch (AGB), intermediate mass stars undergo 
dredge-up and surface enrichment in products of nucleosynthesis. The enriched material is 
thus returned to the ISM via the molecule-rich circumstellar dust shell. 
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More massive (>8 Mq) stars in the red supergiant (RSG) stage may resemble high 
mass-loss rate AGB stars, in that they are also surrounded by optically-thick circumstel- 
lar dust shells, though such high mass-loss rate RSGs are rare. Highly-obscured, evolved 
intermediate-mass (AGB) stars with large mass-loss rates likely dominate the rate of return 
of nuclear processed material to the interstellar medium. However, in our own Galaxy, it can 
be difficult to distinguish between high mass-loss AGB and RSG stars, due to their highly 
uncertain distances. 

IR spectroscopy is particularly useful in determining the properties of circumstcllar 
dust. The emission and absorption properties of different dust grains have been studied in 
the laboratory and modeled (e.g., Draine & Li 2001; Li & Draine 2001). IR-bright stars in the 
Galaxy have been observed with the IRAS LRS and, subsequently, with the spectrometers 
aboard ISO (see reviews by Waters et al. 1999; Waters & Molster 1999; Barlow 1999, 1997 
and references therein) and classified using the dust features revealed in these spectra (e.g., 
Kraemer et al. 2002). The sensitivity of the Spitzer Infrared Spectrograph^ {IRS, Houck et 
al. 2004) on the Spitzer Space Telescope (Werner et al. 2004) now makes it possible to study 
the IR spectra of IR-luminous, evolved stars in nearby galaxies. 

In the absence of spectroscopy, photometric tools such as color-color diagrams can be 
used to identify and classify evolved stars. New studies of stellar populations in the Galaxy 
as well as in nearby external galaxies will rely in particular on photometry obtained with the 
Spitzer Infrared Array Camera {IRAC, Fazio et al. 2004) and Multiband Imaging Photometer 
for Spitzer {MIPS, Rickc et al. 2004). Indeed, the point sources detected in IRAC and MIPS 
imaging of nearby galaxies likely will be dominated by dust-enshrouded post-MS stars. As 
strong silicate and carbonaceous dust features at ~10 //m distinguish C-rich and 0-rich 
envelope chemistries, respectively, the use of IRAC and MIPS colors can, in principal, avoid 
the ambiguity in assigning stellar classes based on near-infrared (NIR) colors (e.g., Egan, 
Van Dyk & Price 2001, hereafter EVPOl). However, such photometric diagnostics must 
first be reliably associated with IR spectral properties, which are best determined from IRS 
spectroscopy. 

As previous infrared surveys have demonstrated (e.g., Groenewegen et al. 1995; Groe- 
ncwegen & Blommaert 1998; Zijlstra et al. 1996; Loup et al. 1997; van Loon et al. 1997, 
1998; Trams et al. 1999), the Large Magellanic Cloud (LMC) is an ideal nearby galaxy to 
study the evolution of mass-losing post-MS stars. Firstly, it contains a large population of 



^This work is based on observations made with the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. The IRS was a 
collaborative venture between Cornell University and Ball Aerospace Corporation funded by NASA through 
the Jet Propulsion Laboratory and Ames Research Center. 
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IR-luminous, mass-losing objects found at essentially the same distance, thereby alleviat- 
ing the distance ambiguities that haunt IRAS and ISO samples of mass-losing stars in the 
Galaxy. This allows the spectral properties of the LMC evolved star population to be directly 
related to their luminosities, providing important constraints on stellar models. Indeed, the 
luminosity distribution of carbon stars in the LMC has long been a key test of models of the 
late stages of stellar evolution in low metallicity environments (e.g., Iben & Renzini 1983). 
Secondly, the low metaUicities and high star formation rates of the LMC mimic those of 
high-redshift galaxies. In addition, spectra of LMC objects can be used to establish "ground 
truth" for IR photometric diagnostics. 

Due to the proximity of the LMC, individual objects have been detected by previous IR 
satellites: IRAS and MSX each detected ~1800 IR sources in or towards the LMC (Schwering 
1989; EVPOl). The most sensitive band on MSX, A band (Ac = 8.3 //m), was approximately 
four times more sensitive than the IRAS 12 /im band, and permitted the detection of evolved 
stars with less extreme IR luminosities than those detected by IRAS (EVPOl). EVPOl 
combined the MSX A band flux (in magnitudes) with near-IR (JHK) measurements from 
the 2MASS survey to form various diagnostic color-color and color-magnitude diagrams for 
the LMC objects. To classify their nature based upon these IR fluxes, they compared these 
color-color and color-magnitude diagrams to analogous diagrams constructed for various 
classes of stars based upon the Wainscoat et al. (1992) "SKY" model of the Galaxy. On 
this basis, they characterized the stars in specific categories, including red supergiants, early 
(low mass loss) 0-rich and C-rich AGB stars (EVPOl category 5), late (high mass loss) 
0-rich and C-rich AGB stars (category 6), and PNs and H II regions (category 7). Within 
these categories, the different types of objects were distinguished somewhat on the basis 
of colors, although only for the late C-rich and 0-rich AGB stars were the objects well 
distinguished within a category. Note that EVPOl classified ~500 of the ~1800 objects 
in their study as foreground Galactic objects. To check the accuracy of this classification 
scheme for the LMC objects, EVPOl investigated the spectral types of the LMC objects 
using data in SIMBAD. However, for the bona fide LMC objects, only a small percentage 
had classifications; these were especially uncommon for categories 6 and 7. Although the 
limited data do not contradict their classifications, these classifications suffer from limited 
spectroscopic support. This is not surprising, as the IR-luminous objects detected by MSX 
do not commend themselves to visible spectroscopy, and the requisite mid-lR spectroscopy 
has been impractical until only recently. To date, the study by EVPOl stands as the most 
comprehensive classification of luminous 8 ^m. sources in the LMC. Therefore it is important 
to verify and/or correct the EVPOl color classifications. 

We have conducted a spectroscopic study with Spitzer IRS (PI: J. Kastner; PID: 3426) 
of a representative sample of objects selected on the basis of their IR colors and magnitudes 
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to be the most IR-luminous post-MS stars in the LMC. Here we present the atlas of spectra. 
We then use these spectra as the basis both for understanding the high-luminosity end of 
the LMC evolved star population and for establishing photometric diagnostics for classifying 
IR stars. In §2 we describe the sample selection, observations, and data reduction. In §3 
we present the spectra and classify them according to continuum and spectral features and 
bolometric luminosities. In §4 we describe the classes of objects identified and compare the 
classifications with the types assigned by EVPOl. We then use our library of IRS spectra 
to show how Spitzer / 2MASS colors can be used to classify IR-luminous sources in external 
galaxies reliably and unambiguously. On the basis of our new classifications, we discuss the 
imphcations for the late stages of stellar evolution. Detailed analysis of individual spectra 
and of the spectral classes established in this study will be presented in future papers. 



2. SAMPLE AND OBSERVATIONS 

The complete sample observed by IRS comprises 62 IR sources, selected on the basis 
of large 8.3 jim fluxes from the comprehensive catalog of 2MASS / MSX sources compiled 
by EVPOl. The EVPOl sources were classified based on IR colors. These sources have 
coordinates accurate to < 1". We selected those MSX 8.3 /im (A-band) sources with Fg.s > 
150 mJy (A< 6.5) that were classified by EVPOl, on the basis of 2MASS / MSX IR colors, 
as "late" (i.e. very red) 0-rich AGB ('O'or 'OH/IR'), "late" C-rich AGB ('C IR'), red 
supergiants ('RSG'), or planetary nebulae ('PN'). This resulted in a preliminary sample 
of ~ 170 sources. Of these, - 50% were classified as OH/IR, - 20% as PN, - 12% as 
RSG, ~ 12% as O AGB, and ~ 5% as C IR. We then randomly deselected sources so as 
to "thin" this preliminary list to 59 targets, distributed so as to include 10 - 15 sources 
from each of these five evolved star spectral classes defined in the EVPOl 2MASS/MSX 
study. To ensure representation in the sample of the most rapidly mass-losing, and hence 
highly obscured, objects — which may have not been detected by 2MASS and therefore may 
have been excluded from the EVPOl sample — we include seven additional IRAS sources. 
These were sources selected from known or candidate evolved stars (Trams et al. 1999), that 
satisfied our 8 /xm fiux criterion, for which 2MASS positions were made available subsequent 
to pubhcation of the EVPOl MSX LMC catalog, and which were not fisted on the Spitzer 
Reserved Observations Catalog. Of this final sample of 66 sources, four targets (1 EVPOl 
source and 3 IRAS sources) were not observed, due to confiicts with other Spitzer programs. 
The 62 observed sample objects are listed in Table 1, where we fist 2MASS (J, H, and K) 
and MSX (A) magnitudes, IR types determined from the IR colors by EVPOl, and object 
types and other star names we obtained from SIMBAD. 
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The targets were observed by Spitzer using the IRS in staring mode. The short-low 
(SL) and long-low (LL) modules were used; these modules cover the wavelength ranges 5.2 
- 14 /im and 14 - 38 /xm, respectively, with resolving powers of 64 - 128. Two nod positions 
along the slit were observed to facilitate rejection of artifacts in the data. For most objects, 
four exposures of 30 and 14 seconds for the SL and LL modules, respectively, were obtained 
at each nod position. A few unusually bright objects were observed with shorter exposure 
times to avoid saturating the detector. To accurately position each object in the sht, a 
peak-up procedure was used {IRS or PCRS). Two observations (MSX LMC 769 and 773) 
were unsuccessful due to peak-up failure; the fields were crowded and the peak-up centered 
on the wrong object. Thus, this atlas presents spectra of 60 objects. 

The raw data were processed through the Spitzer pipeline version Sll.O. After pipeline 
processing, the basic calibrated data (BCDs) were processed using the SMART software^ 
(Higdon et al. 2004). The two-dimensional spectral images (BCDs) of the multiple exposures 
were median-combined and, in most cases, the off-source images were subtracted to remove 
the sky background. In the case of four objects, the off order was contaminated and could 
not be used for sky subtraction. In these cases, the sky was estimated locally, on either side 
of the extracted aperture. 

Spectra were extracted and the fluxes were calibrated using the default point source 
extraction apertures in SMART. The modules at each nod position were merged and the 
edges and overlapping regions of the modules were trimmed. The spectra from the two nod 
positions were then averaged to produce the final spectrum. The Spitzer pipeline produces 
uncertainty images; however the uncertainty images produced by the Sll pipeline are un- 
reliable. Therefore, we estimated the flux uncertainty in each spectrum based upon the 
differences between the two nod position spectra at each wavelength. The flnal uncertainty 
for each spectrum was assumed to be the median value of the uncertainties at all wavelengths; 
these generally ranged from 2-20 mJy, though 4 objects had median uncertainties >30 mJy. 
We found that the relative uncertainties (the median value of the uncertainty at each wave- 
length divided by the flux at that wavelength) were typically ~2%, with a maximum of 5%. 
We note that the wavelength calibration in the Sll pipeline is incorrect (too red) by 0.04 ^m 
for the SL module and also uncertain (by <0.14 /xm) for the LL module. These wavelength 
inaccuracies are not expected to have any effect on the analysis or conclusions of this paper. 



^ SMART was developed by the IRS Team at Cornell University and is available through the Spitzer 
Science Center at Caltech. 
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3. RESULTS: THE SPECTRA 

3.1. SPECTRAL FEATURES 

Mass-losing evolved stars are expected to show distinctive IR features, depending on the 
chemistry of the circumstellar dust (e.g., Draine & Lee 1984). Oxygen-rich dust shows sihcate 
features at 9.7 and 18 /im which appear in emission at low optical depths and absorption 
at high optical depths (i.e. high mass- loss rates). Carbon-rich stars show a broad SiC dust 
emission feature at 11.5 /im and a narrow absorption feature at 13.7 /xm, due to C2H2 gas 
in the circumstellar envelope. 

The IRS spectra are presented in Figure 1. The spectra are grouped according to the 
dominant spectral features: features indicative of a predominance of silicate dust (Fig. la); 
features indicative of a predominance of carbonaceous dust and/or molecules (Fig. lb); 
or polycyclic aromatic hydrocarbon (PAH) features and narrow emission lines (Fig. Ic). 
The 0-rich objects are further divided into subgroups (see §3.3). Within each group, the 
objects are then placed in order of increasing luminosity, which we have calculated (see §3.2). 
Figure 1 also shows the infrared spectral energy distributions (SEDs) of the sample objects, 
including the IRS spectra, MSX 8.3 iim fluxes, and 2MASS and IRAS fluxes (obtained from 
SIMBAD) where available. Two objects (MSX LMC 1072 and 1524) were observed using 
only the LL modules, as data were already being obtained for Spitzer project FID: 1094. A 
third object (MSX LMC 1280) was observed using only the SL modules, as data were being 
obtained for program PID:3505. 

The largest of the three groups of spectra is that of the oxygen-rich objects. These mostly 
show the expected strong silicate emission features and have continua that are dominated by 
photospheric emission at short wavelengths and dust emission at longer wavelengths. The 
silicate features are generally quite smooth, indicating an amorphous composition for the 
silicate dust (MSX LMC 1326 is a notable exception, showing peaks indicative of crystalline 
silicate dust; Kastner et al. 2006). Many of the spectra of the evolved 0-rich objects (most 
of which are RSGs; see §4.1) also show a narrow 6.2 /im emission feature. This is unex- 
pected, since this feature is generally attributed to PAHs, carbon-based molecules which are 
not expected to be found in the circumstellar environment of an 0-rich evolved star. The 
2-dimensionaI spectra show that this feature is commonly seen in emission in the sky (i.e., 
in the LMC) in the vicinity of these 0-rich objects, and that there is some spectral contami- 
nation that is not completely canceled out in the sky subtraction process. However, this sky 
emission is at very low flux density levels, not enough to explain the bright features observed. 
Extraction of off-source spectra indicates that the 6.2 /im feature can is only detectable on 
the source, and is not apparent in off-source spectra, with or without sky subtraction. The 
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presence of this 6.2 /xm carbon feature, and the suggestion of the 11.3 //m PAH feature, 
in the spectra of some 0-rich objects is a puzzle. The second group of objects, the C-rich 
sources, show the expected SiC emission and C2H2 absorption features and have continua 
dominated by dust emission. The broad "30 micron" feature seen in C-rich evolved sources 
and attributed to MgS is also present in most. There are also a significant number of very 
red objects that show forbidden emission lines and strong features attributed to PAHs; these 
make up the third group of sources in our sample. In Figure 2 we present a typical spectrum 
from each of the three groups, with the primary spectral features indicated. 

Spurious features appear in several of the spectra, at wavelengths 15.6, 18.7, 33.5, 34.8, 
and 36.0 /^m. Examination of the 2-dimensional spectral images shows that these are due to 
background "sky" emission lines originating in diffuse gas in the LMC; the lines appear in 
off-source as well as on-source spectral images and are unrelated to the program sources. The 
wavelengths correspond to known [S III], [Si 11], and [Ne 111] emission hnes. These spurious 
features appear in our spectra of C-rich and 0-rich sources in emission or absorption; because 
this diffuse emission is not uniform across the region covered by the spectrograph slit, it is not 
always removed properly by our sky subtraction. These features are seen at much stronger 
levels in emission in the very red emission-line objects, and in these cases they clearly arise 
from the sources themselves. 



3.2. LUMINOSITIES 

Infrared luminosities were calculated by integrating the fiux over the wavelength range 
1 - 100 fim. The IRS spectra were used to determine the fiux densities over the wavelength 
range 5.0 - 36.0 fim. The integrated fiuxes over the range 1-5 fira were determined using 
a blackbody scaled to match the 2MASS fiux densities, and the fiux densities from 36 - 
100 /xm were estimated using a blackbody weighted by A"^'^ (to approximate optically-thin 
emission) and scaled to match the IRS spectrum and available IRAS data. An example of 
this procedure for each of the three groups is shown in Figure 3. A distance to the LMC of 
50.1 kpc (distance modulus 18.50 mag; Freedman et al. 2001) was assumed. 

The uncertainties in the luminosities are dominated by the uncertainty in the shapes 
of the spectral energy distribution in the wavelength ranges outside the IRS spectrum. The 
temperatures of the blackbodies were varied to estimate the uncertainties in the measurement 
of the luminosities, which were found to be 5 - 10% for the 0-rich and C-rich evolved stars. 
The assumption that the shape of the SED between 2 and 5 ^m. can be approximated 
by a blackbody curve was investigated using ISO spectra of Galactic objects and might 
contribute an additional uncertainty to the luminosities of up to a further 19% (typically 



-9- 



4%). In calculating the luminosities of the mass- losing stars, we note that we have neglected 
the effects of variability. Long period variables may vary by up to 1 magnitude in the 
NIR, which results in up to a factor 2.5 uncertainty in the derived luminosities of these 
objects, assuming that the shape of the 1 - 100 /xm SED remains the same as the flux varies. 
Although not affected by variability, the uncertainties in the measurement of the luminosities 
of the red emission-line objects are large because most of the emission is from wavelengths 
longer than the IRS spectral range, where the shape of the spectral energy distribution is 
poorly constrained. The uncertainties in the luminosities of these sources may be up to 
a factor of a few. The bolometric luminosity will be dominated by the IR luminosity for 
most sources. The RSGs are the only stars in the sample for which optical emission might be 
expected to contribute significantly to the SED; however, we find that the derived bolometric 
correction for the RSGs has very little dispersion between objects (§4.4), so there can be 
little contribution shortward of JHK. Luminosities are listed in Table 2. Histograms of the 
luminosities of the stars in each group are shown in Figure 4. 



3.3. CLASSIFICATION OF THE SPECTRA 

Several IR spectral classification schemes have been established since the advent of satel- 
lite mid-lR spectroscopy. IRAS had a Low-Resolution Spectrometer (LRS) that obtained 
spectra in the range 8—23 with a resolution that varied from ~40 at the short end to ~20 
at the long end. The original IRAS LRS atlas contained spectra of 5425 sources (Olnon 
et al. 1986). These were classified based upon the continuum shape (blue or red) and the 
presence or absence of certain dominant spectral features (i.e., silicate emission or absorption 
or C-rich features). Volk & Cohen (1989), Volk et al. (1991), and Kwok, Volk, & Bidelman 
(1997) extracted spectra of additional objects from the IRAS database to arrive at a total 
of 11,224 sources. These were classified into nine groups, with a one-letter designation based 
on the dominant spectral feature(s) and/or the shape of the continuum. 

Kraemer et al. (2002) developed a new classification system (hereafter the KSPW sys- 
tem) that they applied to the over 900 objects for which complete spectra were obtained 

by ISO with the Short Wavelength Spectrometer (SWS). The SWS had a wavelength range 
of 2.4—45 fim with a resolution of 1000 - 2000. Since our IRS spectra extend to 38 fim 
and have a resolution of 64 - 120, they are better matched to the ISO spectra, and we will 
consequently classify our spectra according to the KSPW system. 

We compared our IRS spectra to typical ISO spectra in the KSPW spectral classes^ . In 



^ISO spectra of template objects were obtained from http://isc.astro.cornell.edu/ sloan/library/swsatlas/aotl.html. 
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the KSPW system, stars are classified according to their continua, with subclasses assigned 
according to the appearance of specific spectral features. The continuum is characterized by 
the temperature of the dominant emitter, in particular whether it is photospheric or dust 
emission. Spectral features include broad features associated with silicate or carbon rich 
dust emission, complex features thought to be due to PAH emission, and narrow atomic and 
molecular features. We classified each of our sample objects according to the KSPW system. 

We also assigned an "/i?^ type" to each star in our sample based on the observed spectral 
features and luminosities (as well as optical information for the two peculiar objects). The 
various classifications of the sample objects so derived from the IRS spectra arc listed in 
Table 2, along with the original EVPOl type. The new classifications are discussed in the 
next section. 

4. DISCUSSION 

4.1. O-RICH STARS 

0-rich objects represent the largest group in the sample, with 33/60 objects (55%) 
having IR features indicative of 0-rich dust. The majority of these 0-rich spectra (25/33) 
are dominated by photospheric emission at shorter wavelengths and by dust emission at 
longer wavelengths (group 2 in the KSPW system). Most of these (22 objects) have strong 
silicate dust emission features, and hence we classify these objects as belonging to KSPW 
subgroup 2SEc, while the remaining three arc best classified as 2SEa. These 3 objects show 
weaker, broader 10 fira emission and weaker, more sharply peaked 18-20 /xm emission than 
the 2SEc stars, as well as a 13 /xm feature in emission, possibly due to aluminum-rich silicate 
grains (Sloan et al. 2003). All except one (MSX LMC 1677, discussed below) of the 22 0-rich 
objects with strong dust features and significant photospheric emission have IR luminosities 
> 5 X 10^ L0, and so we identify them as RSGs. The previous classifications of these objects 
(EVPOl) were either 0-rich AGB stars or RSGs (see Table 2). 

Examination of the luminosities of the three 0-rich stars in subgroup 2SEa and MSX LMC 1677 
(subgroup 2SEc) reveals that, if these stars were at the distance of the LMC, they would 
have luminosities >5.8xl0^ L©. These luminosities would place them significantly above the 
distribution of luminosities for the RSGs (Figure 4; these four objects are not included in 
the histogram). We conclude that these four stars are in fact in the halo of the Milky Way 
Galaxy. One of the stars, MSX LMC 412 (RS Men), is a known Galactic Cap Mira variable 
with a heliocentric radial velocity of >140 kms~^ and period of 304 days (Whitclock et al. 
1994). These authors derive a distance to this star of 4.75 kpc using the period-luminosity 
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relation to determine the luminosity and inferring the distance from the measured flux. We 
estimate luminosities for the other three Galactic stars using their K magnitudes and color- 
period- hmiinosity relationships determined by Whitelock et al. (1994). Distances of 4.2, 3.3, 
and 5.1 kpc, for MSX LMC 1150, 1686, and 1677, respectively, were inferred for the three 
sources using luminosities so derived. Three of the stars are identified by EVPOl as O AGB 
stars from their blue NIR colors, while one star, MSX LMC 1686, had been misidentified 
as a C AGB star. The NIR colors of MSX LMC 1686 are consistent with the categories 
of both O AGB and C AGB (EVPOl) and EVPOl note that there is some confusion and 
overlap between these categories. These Galactic Cap stars will be discussed in detail in a 
forthcoming paper. 

Three, possibly four, of the 0-rich stars, MSX LMC 642, 1072, 1524, and 1280, have 
spectra dominated by warm dust (group 3). We classify these as KSPW subgroup 3SE, i.e., 
similar to the majority of the 0-rich stars (subgroup 2SEc) with strong silicate emission at 10 
and 18 /xm from optically thin diist emission, but with thermal continuum from an optically 
thick dust shell. The SED, F^,, peaks around 10 - 15 jim, indicating a strong warm dust 
contribution, although vFy peaks closer to 2 /xm, due to the contribution from the reddened 
stellar photosphere. The silicate emission features are stronger than the Galactic 0-rich 
AGB stars discussed above. The combined SL-I-LL spectra of MSX LMC 1072 and 1524, 
kindly provided by Marckwick- Kemper, confirm the classifications we made based on our LL 
spectra. The classification of MSX LMC 1280 is based on our SL spectrum, which covers 
only the 5-14 //m wavelength range, and is therefore somewhat tentative. The luminosity 
of this source was derived by assuming the same spectral shape as MSX LMC 1189, and is 
highly uncertain. These four 0-rich stars with dust-dominated spectra have IR luminosities 
of ~2 - 7 xlO"^ L0, placing them very high on the AGB, so we tentatively identify these as 
luminous 0-rich AGB stars. Although the theoretical limit for the luminosity of AGB stars is 
6 X 10^ L0, luminosities slightly higher than this may be possible from higher mass (~10 Mq) 
progenitors at low mctallicity (i.e., super- AGB stars Herwig 2005). Our classifications of 
these objects based upon the mid-IR spectra differ significantly from those of EVPOl: One 
of the 3SE stars (MSX LMC 1072) has EVPOl type PN, but its IRS spectrum clearly lacks 
the red continuum. It is likely to be a O AGB star, as it only narrowly falls outside of the 
criteria for identifying O AGB stars in EVPOl: H-K = 0.76, 0.01 magnitudes too red. The 
other stars classified as 3SE have EVPOl type C IR. MSX LMC 1280, 1524, and 642 have 
2MASS/MSX colors inconsistent with the EVPOl criteria for O AGB, OH/IR stars and 
RSGs but their IRS spectra clearly indicate 0-rich envelopes. 

We find no population of less luminous 0-rich AGB stars in the sample. It is possible 
lower luminosity 0-rich objects were excluded from the sample due to our A-band selection 
criterion (Fg.s > 150 mJy; §2). Unlike the C-rich AGB stars, the 8 iiva. fiux is a small 
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fraction of the bolometric luminosity for the 0-rich AGB stars, therefore lower luminosity 
0-rich AGB stars are relatively fainter at 8 fim. The 4 Galactic 0-rich AGB stars in our 
sample, for example, would have Fg.a < 40 mJy if they were at the distance of the LMC, and 
so would have been excluded from our sample. Identifications of 0-rich AGB stars fainter 
than those studied in this work are necessary to confirm the possibility that the lack of lower 
luminosity 0-rich AGB stars in our sample is a selection effect. 

Two of the 0-rich stars show self-absorbed silicate emission features at 10 fim, indicative 
of optically thick dust shells. One of the objects, MSX LMC 1182, shows broad silicate 
emission with strong, narrower central absorption, and has an IR luminosity of 4.7x10^ L©, 
confirming its status as an OH/IR supergiant (Wood et al. 1986; Roche et al. 1993). The 
spectrum of MSX LMC 1182 is dominated by very cool dust so we classify it as group 5 in 
the KSPW system, though its exact subclassification is ambiguous. This star is one of the 
IRAS sources included in our sample to ensure representation of the most highly obscured 
objects (§2). It was not selected on the basis of a classification by EVPOl; it is in the EVPOl 
sample, but does not have an EVPOl type. The other object, MSX LMC 807, shows less 
strongly self-absorbed 10 fim emission, and has luminosity 4x10'^ L0. We identify this object 
as a likely OH/IR star on the AGB, in agreement with its EVPOl type, although we note 
that OH maser detection is required to confirm its formal OH/IR star classification. This 
spectrum falls into KSPW subclass 3SB, which is associated with AGB and OH/IR stars. 

The remaining two 0-rich objects, MSX LMC 890 and MSX LMC 1326, show relatively 
fiat continua with strong broad silicate emission at 10 fim and 18-20 fim. We assign these 
objects IRS type "Pecuhar". MSX LMC 890 and MSX LMC 1326 are optically identified 
as hypcrgiant B stars R 126 and R 66, respectively, and the IR spectra are best matched 
by the class 5SE template spectrum, which are generally young stars. The IRS spectra 
of these stars are likely to be indicative of circumstellar, perhaps circumbinary, disks, and 
are modeled and discussed in detail by Kastner et al. (2006), who classify them as B[e] 
hypergiants. MSX LMC 890 and MSX LMC 1326 have EVPOl type PN based on their NIR 
colors. EVPOl separate PNe from early type stars with surrounding dust (such as these B 
stars) using the H-K color, as PNe are redder than the younger objects in the Wainscoat et 
al. (1992) model. They note that their classification may be confused for objects with H-K 
colors near the boundary, i.e. H-K = 0.75, and other authors have shown that Galactic PNe 
can have NIR colors considerably bluer than this (Whitelock 1985; Ramos-Larios & PhiUips 
2005). MSX LMC 890 and MSX LMC 1326 are slightly redward of this boundary (H-K 
= 0.94 and 0.81, respectively), and so were misclassified by EVPOl as PNe when they are 
instead dusty B stars. It has been noted by Whitelock (1985) that the H-K colors of PNe 
are in fact not dissimilar from hot stars surrounded by dust, though a J-band excess in PNe 
due to He line emission produces redder J-K colors which may distinguish these objects from 
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younger stars. 

4.2. C-RICH STARS 

Sixteen of the 60 sample stars (27%) have spectra typical of warm carbon-rich dust 
envelopes. The IR luminosities of the C-rich stars are 5x10^ - 2x10^ Lq, so we identify all 
16 C-rich stars as AGB stars. The large number of C-rich AGB stars in the sample stands in 
stark contrast to the lack of 0-rich AGB stars found (§4.1) and has imphcations for the late 

stages of stellar evolution (sec §4.5). In addition to the silicon carbide emission feature and 
the narrow C2H2 absorption feature at 13.7 /xm, there is a C2H2+ HCN absorption feature at 
7.5 /xm seen in almost all these spectra. Another broad absorption feature, due to C2H2 and 
HCN, may be present between 14 - 16 fim (Sloan et al. 1998; Aoki et al. 1999), particularly 
in MSX LMC 1384. A broad emission bump, peaking at wavelengths ~26 /^m, is present in 
11 objects. This feature, known as the "30 micron" emission feature, is usually attributed 
to MgS (Goebel & Moseley 1985). Forrest et al. (1981) detected it initially in the most 
optically thick Galactic carbon stars and in C-rich PNe. It has subsequently been seen in 
several carbon-rich PNe and approximately a dozen carbon-rich proto-PNe (Volk et al. 2002). 
There is evidence from ISO that this broad feature is composed of two features, peaking at 
~26 ^m and ~33 fim, and with varying strengths depending upon the evolutionary state of 
the star (Hrivnak et al. 2000; Volk et al. 2002). Hony et al. (2002) have examined this broad 
feature in a variety of evolved carbon-rich objects and are able to fit it by using MgS grains 
of varying temperatures. 

Most of the carbon-rich stars (12/16) have spectra that are dominated by warm dust. 
These are classified as subgroup 3CE. Only three of these 12 have EVPOl type C IR, while 
one has no EVPOl type (no 2MASS data available) and 8 have EVPOl type OH/IR. These 
eight stars have very red NIR colors for carbon-rich stars, compared with the Wainscoat 
et al. (1992) model predictions and EVPOl classifications. Only one star (MSX LMC 218) 
has H-K and K-A colors close to the boundary for C IR stars, while the other seven have 
colors too red to be categorized as EVPOl type C IR, C AGB or RSG. Two carbon-rich 
stars, MSX LMC 1652 and 1592, have spectra dominated by cooler dust; they fall into 
subgroup 3CR. MSX LMC 1652 has EVPOl type OH/IR, with extremely red NIR colors. 
MSX LMC 1592 has no 2MASS data so has no EVPOl type. One star, MSX LMC 1400, 
shows clear SiC features, but matches equally well the 3CE and 3CR templates, so has been 
classified as '3C?'. It also has very red NIR colors for a C-rich star and has EVPOl type 
OH/IR. 

The final object classified as a carbon-rich star, MSX LMC 775, shows a nearly feature- 
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less spectrum with a narrow peak at 5 - 6 //m and relatively blue continuum. The 13.7 jim 

C2H2 absorption feature apparent in all other carbon-rich stars in our sample is notably 
absent. This spectrum is not well matched by any of the naked star (group 1) or stellar 
photosphere + dust (group 2) ISO template spectra. The spectrum of MCX LMC 775 is 
very similar in slope to that of MSX LMC 1492, but without the clear 13.7 feature. The 
lack of strong 0-rich or C-rich dust features suggests that this object may be an S star, with 
approximately equal carbon and oxygen abundances (e.g.. Jura 1988). However, the IRS 
spectrum of MSX LMC 775 does not resemble the ISO spectrum of the S star W Aql and 
shows a hint of a SiC feature at 11.5 //m, which is rare for S stars (Chen & Kwok 1993), sug- 
gesting the emission may be due to optically thick carbon-rich dust. This unusual object may 
be similar to R Corona Borealis, which typifies a class of stars that have hydrogen- deficient 
and carbon-rich atmospheres and are highly variable due to periods of dust formation (see 
review by Clayton 1996). These objects show nearly featureless mid-lR spectra, that are 
dominated by blackbody continuum (Clayton 1996; Kracmcr ct al. 2005). Two such R CrB 
stars have recently been identified in the Small Magellanic Cloud (Kraemer et al. 2005); 
these objects show similar Spitzer spectra to MSX LMC 775, although MSX LMC 775 is 
well-described by a modified blackbody (multiplied by A~^), rather than a simple blackbody 
continuum as seen in these two. RCB stars are identified mainly by the absence of hydrogen 
features in their optical spectra and by variability characterized by dramatic declines in the 
optical brightness at irregular intervals (Clayton 1996). Further data, such as optical spec- 
troscopy and an optical light curve, are necessary to confirm if this candidate is indeed an 
R CrB source or an S star. MSX LMC 775 has EVPOl type C JR. 



4.3. RED, EMISSION-LINE OBJECTS 

Approximately one fifth of the sample (11 objects) is comprised of very red objects 
which have a peak in their SED longer than the wavelengths observed by IRS (group 5). 
These objects arc dominated by cool dust and show strong PAH emission features and narrow 
emission lines. These red, emission-line objects were classified by EVPOl as type PN, which 
they distinguish from H II regions associated with massive star formation by H-K colors 
(H-K > 0.75). Indeed, the spectra of these objects resemble the ISO spectra of PNe that are 
classified as 4PU in the KSPW scheme. However, various fines of evidence, discussed below, 
indicate that these objects are H II regions rather than PNe. 

The luminosities of the red objects are rather high, 2 x 10*^ - 5 x 10^ L0. Most of these are 
too luminous to be the immediate descendents of evolved, intermediate-mass stars (Figure 
4) , but rather they agree with the luminosities of H II regions powered by massive young O 
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stars. The IR luminosities of many of these objects are somewhat uncertain, as much of the 
flux is emitted beyond 40 ^m, where the shape is the spectral energy distribution is poorly, 
if at all, constrained (§3.2). The infrared flux densities of the H II regions are, if anything, 
more likely to be underestimated, due to extended emission falling out of the slit. Thus the 
luminosities argue against most of these objects being PNe. 

Examination of the optical (Digital Sky Survey) and IR {MSX A band) images reveals 
widespread nebulosity around these objects, suggesting they are the bright, compact cores of 
H II regions associated with massive star formation regions, rather than PNe. This is shown 
in Figure 5. The MSX 8.3 fim. images in particular probe the dust and reveal extended 
nebulous emission, probably dominated by PAH's. Several of the DSS images also show 
diffuse emission, possibly dominated by Ha, while, in contrast, the 2MASS NIR images, 
which we do not display, generally show only the bright point sources. 

The spectra of the red, PAH-dominated objects match KSPW subgroup SUE (Fig. 6), 
and show strong atomic fine-structure fines and PAH features. Galactic objects in this class 
are mostly young objects. In addition, these IRS spectra are lacking many of the high- 
ionization narrow emission lines common in Galactic PNe (Kraemer et al. 2002) and LMC 
PNe (Bernard-Salas et al. 2004), such as [Ne V] 24.3 pim and [O IV] 25.9 /xm, as well as 
[Ne VI] 7.65 //m. This is shown in Figure 6. The IRS spectra of the PAH-dominated objects 
show strong [S III] emission peaks at 18.7 and 33.5 fim and [Ne II] 12.8 ^m, emission lines 
found in both PNe and H II regions (Givcon ct al. 2002). The similarity of the continuum 
shapes of the PN and H II region in Figure 6 suggests that these types of objects will have 
similar mid-lR (5 - 35 /xm) colors and so will be difficult to distinguish photometrically in 
this wavelength range. 

The IRS spectra provide further evidence that the PAH-dominated sources are extended 
H II regions and not PNe. The discontinuity seen at '^14 /im between the spectra formed 
by the short and long wavelength modules (see Fig. 2c) is attributed to the extended nature 
of these sources (a similar phenomenon also affects certain ISO spectra of extended sources; 
see, e.g., Lamers et al. 1996). This "jump" indicates that the spatial extents of the PAH- 
dominated sources are greater than the width of the SL sfit (3.6"), which is smaller than 
the LL slit (10.5"). A spatial extent of 5" at 50.1 kpc corresponds to a physical size of 
25000 AU (1.2 pc), confirming that these sources are unlikely to be PNe. The red, emission- 
line sources also appear extended in 2-D spectral images and this is confirmed by fitting 
Gaussian functions along the spectrum of each source and comparing the FWHM to a similar 
fit to one of the 0-rich point sources (MSX LMC 1150). 

In summary, we find that all of the program objects selected as PNe based on their 
NIR/ MSX colors are, in fact, dusty H II regions associated with luminous young 0/B stars. 
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Thus the color distinction between these two groups made by EVPOl is not supported. This 

in fact confirms their caveat that many of their objects classified as PN may be H II regions 
with enhanced reddening. Revised criteria will need to be developed to distinguish PN from 
H II regions photometrically, and these will probably need to rely upon wavelengths other 
than just the near- and mid-IR. 

4.4. INFRARED COLOR-COLOR AND COLOR-MAGNITUDE 
CLASSIFICATIONS OF LMC SOURCES 

Figure 7 shows IR color-color diagrams for the stars in this study, based upon the 2MASS 
and MSX photometry. The stars of most of the various IRS spectral types separate quite 
cleanly into different regions of this diagram, confirming that IR color criteria can be used to 
classify the spectral types, although this may not be the case for distinguishing PNe and H II 
regions, and PNe may occupy the region between the AGB stars and the H II regions and 
make the separation less clean. As previously noted, however, our IRS spectral classifications 
frequently differ from the original EVPOl classes assigned to these objects. EVPOl indeed 
mention that the 0-rich and C-rich AGB star classifications are rather uncertain due to the 
similar regions in color space that these objects occupy. Other authors have also shown, both 
observationally (van Loon et al. 1998; Zijlstra et al. 2006) and theoretically (Groenewegen 
2006), the difficulty of distinguishing 0-rich and C-rich evolved stars based on their IR 
colors. We nevertheless propose modified IR color criteria, based on the locations of objects 
in Figure 7, subject to the uncertainty due to our sample incompleteness. The isolated 
locations of the two OH/IR stars in our sample in the color-color diagrams suggests that 
other EVPOl objects with 2MASS / MSX colors K-A>4 and J-K>2, but outside the C-rich 
AGB area, are candidate OH/IR stars. The four objects we classify as 0-rich AGB stars 
lie fairly close to each other on the color-color and color-magnitude diagrams, indicating 
that objects with K-A = 1-3 and J-K = 2-3 are candidate 0-rich AGB stars in the 
LMC. Further observations are required to probe the borderline regions where classifications 
overlap or are uncertain and to identify where PNe in the LMC lie on these diagrams. 

Color-magnitude diagrams are also extremely useful in discriminating between spectral 
types of stars (Figure 8). The Galactic and LMC objects separate clearly, with Galactic 
objects having brighter apparent magnitudes due to their relatively modest distances. The 
AGB star luminosity sequence is also evident, with the C-rich objects in the majority and 
the most IR-luminous AGB stars being 0-rich. The most IR-luminous LMC objects in our 
sample are the RSGs and the lone OH/IR supergiant (MSX LMC 1182), the latter of which 
is redder due to the thick dust shell presumably surrounding this source. The lone OH/IR 
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AGB star (MSX LMC 807) has a lower IR luminosity, comparable to the C-rich AGB stars. 
The H II regions, which have J-band magnitudes similar to, but K magnitudes fainter than, 
the C-rich AGB stars (due to their bluer J-K colors; Fig. 7), display the reddest K-A colors 
(Fig. 7). 

One of the primary goals of this program was to determine reliable relations between 
Spitzer IRAC and MIPS photometry and IRS spectral types for compact, infrared-luminous 
objects, especially evolved stars. This will enable the classification and study of evolved stars 
in other nearby galaxies using Spitzer photometry. Therefore, we also derived IRAC and 
MIPS fiux densities for the objects using the IRS spectra and the spectral response functions 
of the IRAC 5.8 and 8.0 /^m bands and the MIPS 24 /xm band. Figure 9 shows the simulated 
IRAC / MIPS colors of our stars as derived from their IRS spectra. Also plotted on Figure 9 
are Galactic stars in each spectral type, as well as PNe, for comparison. The three Galactic 
PNe plotted overlap with the H II regions and appear to fill the gap between the evolved stars 
and H II regions. Clearly, observations of LMC PNe are required to determine whether the 
proposed photometric classifications of evolved stars and H II regions are "contaminated" 
by PNe. The comparison objects were selected at random from the database of Kraemer 
et al. (2002) — subject to the constraint that the ISO spectrum had reasonable signal-to- 
noise ratio — and the synthetic IRAC and MIPS colors were derived from the ISO spectra. 
The Galactic objects fall in the same regions of the diagram as the LMC stars, indicating 
that the IRAC / MIPS color-color diagrams are not very sensitive to metallicity. Therefore 
these diagnostics can be applied to external galaxies without being adversely affected by 
metaUicity differences, at least in the range ~0.5 - 1 Zq. It is possible to construct a large 
number of color-color diagrams using the observed 2MASS / MSX photometry and synthetic 
Spitzer IRAC / MIPS colors. In Figure 10 we present two such color-color diagrams that 
appear to show the best separation of the object types. Based on these two diagrams, we 
propose color criteria to classify luminous 8 fim sources in external galaxies (shown as boxes 
in Fig. 10 and listed in Table 3). Figure 11 shows color-magnitude diagrams derived solely 
from the synthetic Spitzer flux densities of our sample objects. The flgure indicates that 
such IRAC /MIPS color-magnitude diagrams represent an effective means to distinguish 
between evolved stars and H II regions in external galaxies, and should easily identify the 
most luminous, rapidly mass- losing stars (i.e., OH/IR supergiants) in such galaxies. 

IR colors may also be useful to infer the luminosity of evolved stars. Using the measured 
luminosities of the C-rich AGB stars and the RSGs, for which we have a sufficiently large 
number of objects, we investigated the relation between the luminosity and the NIR (K-band) 
magnitude. We defined a quasi- "bolometric correction" to be the ratio of the bolometric 
luminosity to the K-band fiux density. For the C-rich AGB stars, this ratio shows a clear 
trend with the K-A color (Figure 12). This trend can be fitted with an exponential function. 
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of the form Lioi/Fk — a(K — A)** + c, where L},oi is the bolometric (IR) luminosity in units 

of Lgy^ni is the 2MASS K-band flux density in Jy, and (K-A) is the 2-8 /xm color 
in magnitudes, derived from the 2MASS and MSX K and A magnitudes. The best fit 
parameters are a = 0.0015, b = 6.26, and c = 14.9. Using this fit to the points shown in 
Figure 12, the bolometric correction for a source of a given K-A color can be estimated and 
hence the luminosity derived from the K-band magnitude. The largest error in the luminosity 
predicted in this way is less than a factor of 2, with most of predicted luminosities within 
20% of the measured luminosity. 

For the RSGs, there is a unique value of the "bolometric correction" for all K magnitudes 
(Figure 13). The mean value for all sources is L^oi/Fk = 1.46 ^o.i7 x 10^ Lq/Jj. Thus the 
luminosities of RSGs can be predicted from their K magnitudes with an uncertainty of less 
than 15%. 

These results are potentially very useful in the study of Galactic carbon stars and 
RSGs, assuming that metallicity does not affect the bolometric correction. The bolometric 
corrections derived here will enable estimates of luminosities, and hence distances, of stars 
to be made simply from their IR colors and fluxes. 

4.5. IMPLICATIONS FOR LATE STAGES OF STELLAR EVOLUTION 

The predominance of C-rich stellar envelopes among the lower-luminosity members of 
our sample, and a relatively sharp dividing line between C-rich and 0-rich infrared-luminous, 
evolved stars, is apparent in Table 2 and Fig. 4. These results suggest constraints on models of 
the processes of nuclear burning and dredge-up in intermediate- to high-mass stars and, at the 
same time, may indicate the importance of the metallicity of the star-forming environment 
on the ultimate evolutionary endpoints of such stars. 

We caution that our IR-luminous LMC sample was selected based on a classification 
scheme that has proven somewhat unreliable, especially with regard to the distinctions be- 
tween 0-rich and C-rich AGB stars and between PNe and H II regions. This sample therefore 
may not well represent the IR-luminous, evolved stellar population of the LMC. The flux 
limit imposed may also have biased the sample against 0-rich AGB stars. Nevertheless, the 
fact that we find ~ 80% of the AGB stars in our sample are C-rich argues that the low 
metallicity of the LMC serves to facilitiate carbon star formation, due to the relative ease 
of inverting the C/0 ratio at the stellar surface. At the same time, the apparent upper 
cutoff in the luminosity of IR-luminous carbon stars in the LMC, ~ 2 x 10^ Lq (Table 2), 
implies an upper limit of ~3.7 Mq for the initial mass of carbon star progenitors in the LMC 
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(Marigo et al. 1999), though we note that this depends on assumptions about the mass loss, 
dredge-up efficiency and envelope burning. The apparent presence of such a carbon star 
luminosity cutoff, and the fact that the only three AGB stars brighter than this are 0-rich, 
suggests that "hot bottom burning" (see Herwig 2005 and references therein) inhibits carbon 
star formation in the LMC for initial masses larger than ~3.7 Mq. Furthermore, the fact 
that all RSGs in the sample are 0-rich suggests that the He-burning phase for such massive 
(> 8 Mq) stars is short, that the onset of C burning in such stars is rapid, and/or that the 
envelopes of such stars are sufficiently massive that it is difficult to invert the C/0 ratio even 
given low initial metallicity. 

We also note that, for the 0-rich evolved stars, all except one (MSX LMC 1072) with 
LiR> 6 X 10^ L0 have spectral shape 2SEc, while all except one (MSX LMC 815) with 
LiR< 6 X 10^ Lq have spectral class 3SE. This transition coincides with the theoretical upper 
limit on AGB luminosity (Iben & Renzini 1983), although luminous 0-rich stars in a hot 
bottom burning phase may fall above the classical AGB limit (e.g.. Blocker & Schoenberner 
1991) so the luminosity separation of AGB stars and RSGs is uncertain. The two exceptions 
are the ones closest to this dividing luminosity and it may be that variability is playing a role 
in confusing the classification of the two objects. Thus is it possible that the IRS spectral 
shape may provide clues to an object's luminosity and progenitor mass, which would be 
particularly useful in cases where the luminosity is otherwise ambiguous. 

We caution, however, that the present LMC IRS sample was selected to represent the 
types of object in the LMC but not their relative numbers, and so is not proportionally 
representative of the luminous mid-IR population of the LMC. For example, the most lu- 
minous sources hsted in EVPOl are dominated by objects classified as "OH/IR stars" — 
which, our IRS survey demonstrates, are predominantly or perhaps exclusively C-rich AGB 
stars — whereas such objects constituted only ~ 20% of the sample we have surveyed thus 
far with IRS. Thus, it is important to revisit the complete sample of luminous 8 /xm sources 
in the EVPOl lists, to reclassify these objects according to our revised criteria and, thereby, 
determine the distribution of envelope spectral types among these stars. This work will be 
the subject of a future paper. In addition, it is crucial to obtain additional IRS spectra 
of the most highly obscured yet luminous evolved stars. Such objects would likely appear 
as MSX and IRAS sources but not as 2MASS sources and, hence, are not necessarily well 
represented in the present sample. The OH/IR supergiant MSG LMC 1182, along with the 
two luminous G-rich AGB stars that are also IRAS sources (MSX LMC 1592 and 1652), 
offer indications of the types of sources we might expect to find among this group. 
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5. CONCLUSIONS 

We have obtained Spitzer low-resolution spectra for a representative sample of 60 lumi- 
nous 8 /im soiuxes in the LMC. The sources observed were selected to include a variety of 
evolved, mass- losing stars, based on recently established 2MASS / MSX color criteria. The 
IRS spectra covered the wavelength range 5-35 /xm, allowing determination of the envelope 
chemistry from detailed spectral features, particularly silicate emission features at 10 and 
18 fj,m in 0-rich objects and the SiC emission feature at 11 /im and C2H2 absorption feature 
at 13.7 //m in C-rich objects. We derived IR luminosities from the combined 2MASS and 
IRS spectral energy distributions of the objects and used these, along with the spectra, to 
classify the object types. We identify 16 C-rich AGB stars, 4 0-rich AGB stars, 21 RSGs, 
11 H II regions, and 2 OH/IR stars (one supergiant and one AGB star). We also find 2 B 
supergiants that have peculiar IR spectra and 4 0-rich AGB stars that are in fact foreground 
Galactic Mira variables. We find the vast majority of the AGB stars, with IR luminosities 
~ 10^ Lq, have C-rich envelopes, while the 0-rich objects are almost all luminous RSGs 
with LiR> 10^. The predominance of C-rich stars along the AGB in the LMC is consistent 
with the hypothesis that C-rich stars form more easily than 0-rich stars in lower metaUicity 
environments. However, all C-rich AGB stars in our sample have Lir < 2 x 10^ L©, while 
all four 0-rich AGB stars have Lir > 2 x 10^ L©, suggesting a relatively firm upper limit 
(~3.7 Mq) for the progenitor mass of LMC Carbon stars. The fact that all of the supergiant 
stars are 0-rich, in turn, may indicate that massive supergiants do not have sufficiently long 
He-burning lifetimes to produce C-rich surfaces, and hence C-rich ejecta, despite the low 
metaUicity of the LMC. 

Comparison of the IRS classifications with the NIR color-based types (EVPOl) revealed 
a number of inaccuracies in the 2MASS / MSX classification criteria. All of the sample objects 
with EVPOl type PN were reclassified, mostly as H II regions. Stars classified as OH/IR 
stars on the basis of NIR colors are mostly C-rich AGB stars, based on their IRS spectra. 
The 0-rich AGB stars in the sample were reclassified as either RSGs or Galactic O AGB 
stars, based on their IRS spectra and luminosities. On the basis of these reclassifications, we 
have developed revised 2MASS / MSX color-color and color-magnitude classification criteria 
for LMC objects and, by extension, for luminous 8 /xm sources in other external galaxies. 

For the C-rich AGB and RSG stars, bolometric corrections to the stellar K-band flux 
densities were derived using the observed IR SEDs. For C-rich AGB stars, the bolometric 

corrections were found to depend on the K-A color, allowing the bolometric luminosity to 
be estimated from the K and A-band fluxes with a typical uncertainty of less than 20%, and 
at most a factor of 2. The bolometric corrections of the RSGs were found to be independent 
of IR color; for these objects the bolometric luminosities can be predicted from the K-band 
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magnitude with an uncertainty of less than 15%. The derived bolometric corrections could 
prove very useful for determining the luminosities of carbon stars and RSGs in the Galaxy 
and other nearby galaxies. 

Infrared {IRAC / MIPS / 2MASS) colors were determined from the IRS spectra and 
2MASS fluxes. We found that the spectral types separate well in various mid-IR color- 
color diagrams, indicating that Spitzer photometry can be an effective means to infer the 
spectral class of dust-enshrouded objects in external galaxies. We suggest diagnostics for 
classifying IR-luminous sources in the Galaxy and nearby galaxies, using Spitzer photome- 
try. These colors may also prove useful in the analysis of Spitzer photometry of unresolved 
stellar populations in high-redshift galaxies. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under a con- 
tract with NASA. Support for this work was provided by NASA through awards issued by 
JPL/Caltech to the authors' institutions. We would hke to thank Dan Watson for providing 
peak-up stars, without which many of these observations would not have been possible. We 
thank the referee, A. Zijlstra, for helpful comments that improved the paper. 
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Table 1. LMC Spitzer sample 



MSX RA" Dec" J H K A EVPOl SIMBAD 2MASS IR.AS Other 

LMC J2000 J2000 (mag) (mag) (mag) (mag) type type*" Name Name Name'' 

No. 



22 


05 


04 


47 


.14 


-66 


40 


30, 


,4 


14.18 


13.75 


12.96 


6 


.06 


PN 


IR source 


J0504471,5-6640307 


05047-6644 




87 


05 


10 


19 


.63 


-69 


49 


51, 


2 


14.60 


12.73 


11.07 


6 


.42 


OH/IR 




J05101962-6949514 






95 


05 


10 


00 


.00 


-69 


56 


09, 


,6 


13.39 


11.49 


10.01 


6 


.27 


C IR 




J05095999-6956097 






141 


05 


05 


33, 


.48 


-70 


33 


46. 


,8 


8.84 


7.97 


7.63 


6, 


.36 


RSG 


M 


J05053350-7033469 




[L72] LH 24- 15 


217 


06 


13 


24, 


.67 


-69 


10 


48. 


,4 


14.69 


14.21 


13.23 


5, 


.92 


PN 


IR source 


J0B13246 5-6910480 


05137-6914 




218 


06 


13 


16, 


.39 


-68 


44 


10. 


,0 


13.39 


11.47 


9.93 


6, 


.15 


OH/m 




J05131640-6844099 






220 


06 


12 


32 


.06 


-69 


15 


40, 


,7 


12.91 


11.15 


9.78 


6 


.21 


C IR 


V* 


J05123206-6915404 




DCMC J051232.11-691B40.B 


222 


06 


13 


42 


.48 


-69 


35 


21, 


,8 


14.74 


15.14 


13.95 


6 


.44 


PN 




J05134244-6935219 






264 


05 


14 


49 


.75 


-67 


27 


19, 


8 


8.62 


7.75 


7.39 


5 


.85 


RSG 


M3Iab: V* 


,105144972-6727197 


05148-6730 


HV 916 


412 


05 


15 


41 


.28 


-73 


47 


13, 


9 


6.59 


5.71 


5.26 


4 


,57 


O AGB 


Mc V* 


,1051,54126-7347137 


05169-7350 


R,S Men 


438 


05 


25 


19, 


.51 


-71 


04 


02. 


,6 


16.06 


13.23 


10.96 


6 


.09 


OH/IR 


IR source 


J05251951-7104027 




LI-LMC 1028 


629 


06 


23 


43, 


.63 


-66 


41 


59. 


,6 


9.00 


8.10 


7.74 


6, 


,24 


RSG 


M3/M4 V* 


J05234361-6641696 


05235-6644 


HV 12793 


649 


06 


26 


11 


.38 


-66 


12 


10. 


,8 


9.32 


8.26 


7.71 


6, 


.04 


O AGB 


Star in cluster 


J05261135-6612111 




NGC 1948 WBT 54 


559 


06 


25 


49, 


.25 


-66 


15 


08. 


,3 


14.38 


13.58 


12.49 


5, 


.98 


PN 


UV source 


J05254923-6615087 




[HCB95] LH 52 4978 


587 


06 


31 


04 


.15 


-69 


19 


03, 


,3 


9.12 


8.12 


7.65 


5 


.64 


O AGB 




J05310418-6919030 






589 


06 


26 


34, 


.80 


-68 


51 


40. 


,0 


8.50 


7.65 


7.29 


5, 


.90 


RSG 


M2Iab:-|-... 


J05263479-6851400 




[GMP94] 301 


593 


06 


28 


28, 


.87 


-68 


07 


08. 


,0 


8.48 


7.66 


7.34 


5, 


.95 


RSG 


MOIa V* 


J05282886-6807078 




HV 2561 


597 


06 


29 


42, 


.24 


-68 


57 


17. 


,3 


8.00 


7.24 


6.96 


6, 


.21 


RSG 


Mlla 


J05294221-6867173 


05300-6869 




642 


06 


28 


48, 


.14 


-71 


02 


28. 


,7 


11.31 


9.70 


8.64 


6, 


.30 


C IR 


M8 


J0B284817-7102289 


05294-7104 




764 


06 


32 


52, 


.68 


-69 


46 


22. 


,8 


13.92 


13.76 


12.84 


4, 


.93 


PN 


IR source 


J0532B272-6946226 


05333-6948 




769 


05 


31 


42 


.41 


-68 


34 


53, 


7 


15.08 


14.55 


13.45 


5 


.91 


PN 




J05314242-6834539 






773 


05 


35 


41 


.14 


-69 


11 


59, 


,6 


9.58 


8.55 


8.06 


6 


.07 


O AGB 




J05354110-6911596 






775 


05 


32 


56, 


.18 


-68 


12 


49. 


,0 


13.24 


11.32 


9.85 


6, 


.44 


C IR 




J05325618-6812487 






807 


06 


32 


37, 


.15 


-67 


06 


56. 


,5 


15.50 


13.63 


11.93 


5, 


.66 


OH/IR 




J0B323716-6706664 






810 


06 


30 


20, 


.66 


-66 


53 


01. 


,7 


8.78 


7.92 


7.72 


6, 


.21 


RSG 


Double star 


J0B302067-6663018 




CCDM J06303-66B3B 


815 


06 


36 


14, 


.09 


-67 


43 


55. 


,6 


9.49 


8.68 


8.14 


6, 


.37 


O AGB 


M4 V* 


J0B351409-6743668 




HV 1001 


836 


06 


32 


31, 


.99 


-66 


27 


15. 


1 


15.94 


16.01 


14.00 


5, 


.85 


PN 


IR source 


J05323195-6627164 


05325-6629 




839 


05 


31 


36 


.82 


-66 


30 


07, 


,9 


8.44 


7.66 


7.36 


6 


.27 


RSG 




,705313681-6630076 






870 


05 


35 


28 


.32 


-66 


56 


02, 


,4 


8.45 


7.65 


7.30 


6 


.20 


RSG 


M3Iab: 


,105352832-6656024 


05354-6657 


Dachs LMC 2-16 


889 


05 


38 


31, 


.63 


-69 


02 


14. 


,6 


14.53 


13.80 


12.68 


4, 


.65 


PN 


Molecular cloud 


J05383167-6902146 




[JGB981 30 Dor-06 


890 


06 


36 


25, 


.87 


-69 


22 


55. 


,9 


10.12 


9.73 


8.79 


4, 


.75 


PN 


B;e Em. object 


J0B362B86-6922668 


05368-6924 


HD 37974 


894 


06 


38 


44, 


.71 


-69 


24 


39. 


,6 


13.91 


14.00 


12.85 


6, 


.26 


PN 


Em. object 


J0B384470-6924396 




LHA 120-N 168B 


897 


06 


40 


43, 


.73 


-69 


21 


57. 


,9 


8.87 


7.96 


7.49 


6, 


.29 


O AGB 




J0B40437B-6921681 






934 


06 


39 


IB, 


.86 


-69 


30 


38. 


,5 


14.84 


14.42 


13.48 


4, 


.95 


PN 


Molecular cloud 


J0B391B87-6930384 




[JGB98] N1B8-2 


939 


05 


40 


48 


.48 


-69 


33 


36, 





9.24 


8.26 


7.75 


5 


.71 


O AGB 


Open cluster 


J0B4048B0-6933360 




[HS66] 386 


943 


05 


41 


10 


.63 


-69 


38 


03, 


8 


8.80 


8.00 


7.66 


6 


.42 


RSG 


MO V* 


J0B41 1066-6938040 






1007'^ 


04 


28 


30, 


.17 


-69 


30 


50. 


,0 


16.29 


13.83 


11.87 


6, 


.66 


OH/IR 






04286-6937 




1072 


04 


40 


28, 


.51 


-69 


55 


13. 


,8 


10.16 


8.92 


8.16 


5, 


.35 


PN 


M7.6 


J04402848-6966136 


04407-7000 






04 


47 


16, 
























IR, source 


JU44 ( lDUy-Do.Z4l.ZOO 






1132 


04 


49 


22, 


.46 


-69 


24 


34. 


,6 


9.11 


8.18 


7.78 


6, 


.17 


O AGB 


M V* 


J04492246-6924344 




HV 2236 


1150 


04 


39 


23, 


.66 


-73 


11 


02. 


,8 


6.92 


5.99 


5.52 


4, 


.70 


O AGB 




J04392369-7311028 






1182" 


04 


55 


10, 


.49 


-68 


20 


29, 


,8 


9.25 


7.74 


6.85 


2 


.43 




M7.5 


.704551048-6820298 


04553-6825 




1189 


04 


55 


03 


.07 


-69 


29 


12, 


8 


8.67 


7.68 


7.20 


5 


.32 


O AGB 


IR. source M2 


,704550307-6929127 


04553-6933 




1204 


04 


55 


16, 


.03 


-69 


19 


12. 


,0 


8.55 


7.74 


7.37 


6, 


.35 


RSG 


Em. -line star M 


J04551604-6919120 




[BB74] 164 


1280 


06 


00 


19, 


.03 


-67 


07 


58. 


1 


12.03 


10.44 


9.28 


6, 


.44 


C IR 


M9 


J0B001899-6707580 


05003-6712 




1282 


06 01 


.86 


-67 


35 


23. 


,6 


17.24 


14.61 


12.18 


6, 


.43 


OH/IR 




J05010087-6736236 






1306 


04 52 


58 


.80 


-68 


02 


56, 


8 


15.85 


15.45 


14.15 


6 


.,34 


PN 


Em. object 


,704525878-6802569 




LHA 120-S 5 


1326 


04 


56 47 


.04 


-69 


50 


24, 


7 


10.05 


9.65 


8.84 


4 


.85 


PN 


Em. -line star B8Ia 


,704564705-6950247 


04571-6954 


HD 268835 


1328 


04 


57 43 


.30 


-70 


08 


50, 


3 


8.42 


7.62 


7.32 


.5 


,60 


RSG 


M4 V* 


J04574331-7008503 


04581-7013 


HV 2255 


1330 


04 


55 


21 


.65 


-69 


47 


16, 


,8 


8.87 


7.95 


7.61 


6 


,03 


RSG 


M 


,70455216,5-6947167 




[M2002] LMC 24410 


1384 


05 


43 36, 


.05 


-70 


10 


35. 


,0 


16.37 


13.53 


11.44 


6, 


,37 


OH/IR 


V* 


J05433602-7010361 






1400 


06 40 20, 


.62 


-66 


14 


44. 


2 


16.60 


13.91 


11.77 


6, 


,02 


OH/IR 




J054020B7-6614442 






1429 


06 


44 


13, 


.70 


-66 


16 


44. 


,8 


9.01 


7.92 


7.48 


6, 


,48 


O AGB 


M0.5 V* 


J05441373-6616446 




HV 2834 


1488 


05 50 6 


.72 


-71 


46 


03, 


,0 


15.14 


12.91 


11.04 


6 


,50 


OH/IR 


IR source 


J05500676-7146026 


05508-7146 




1492 


05 49 8 


.83 


-71 


32 


07, 


1 


12.29 


10.66 


9.49 


6 


,46 


C IR 




J05490888-7132069 






1524 


05 


55 


21 


.05 


-70 


00 


03, 


2 


11.97 


10.32 


9.11 


6 


,04 


C IR 




J05552103- 7000030 


05558-7000 




1592*= 


05 


56 


38 


.76 


-67 


53 


34, 


,4 








5 


,81 




Carbon star 




05568-6753 




1652 


06 02 


31, 


.03 


-67 


12 


46. 


,8 


17.85 


15.64 


12.99 


6, 


,14 


OH/IR 


Carbon star 


J06023105-6712469 


06025-6712 




1677 


06 01 


27, 


.82 


-66 


05 


23. 


,3 


6.81 


5.88 


5.15 


3, 


,05 


O AGB 


IR source 


J06012780-6606231 


06013-6606 




1686 


06 06 47 


.81 


-66 


48 


12, 


,6 


6.33 


6.43 


4.84 


3 


,58 


C AGB 




J06064779-6648126 






1794 


i)r, 


41) 


4.'i 


,<)!) 


-69 




.-)4 




l."),S<) 


i.'').:i6 


14,2.'-! 


(i 


,44 


PN 




J054()44()l-6!)25546 






()4:!74-(iS:ii' 


■''(J4 














:! 


















()4:i74-lj,-;:!i 





^Source positions are from the 2MASS catalog. 
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Source types and names from SIMBAD, simbad.u-strasbg.fr/sim-fid.pl. Spectral types are indicated where known, blank indicates that the object is a star of 
unknown spectral type. V* indicates a variable star, and "Em." indicates emission. 

"^Additional IRAS source included in sample (see §2). 

*^This object was not detected by MSX, so is designated by its IRAS name. IRAS flux densities for this object are F{12 nm): 0.199 Jy; F(25/im): 0.25 Jy. 
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Table 2. Properties of the IRS spectra. 



MSX 


Envelope 


IRS 


LiR 


KSPW 


EVPOl 


Name 


chemistry 


type 


Lq 


class 


type 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


1280^ 


O-rich 


O AGB 


1.7 XIO*: 


3SE 


C IR 


1524 


O-rich 


O AGB 


4.1 xlO"' 


3SE 


C IR 


642 


O-rich 


O AGB 


4.2 xlO* 


3SE 


C IR 


1072 


O-rich 


O AGB 


6.5 XlO* 


3SE 


PN 


815 


O-rich 


RSG 


5.5 xlO" 


2SEc 


O AGB 


1132 


O-rich 


RSG 


7.7 xlO* 


2SEc 


O AGB 


549 


O-rich 


RSG 


7.9 xlO* 


2SEc 


O AGB 


529 


O-rich 


RSG 


7.9 xlO* 


2SEc 


RSG 


939^ 


O-rich 


RSG 


8.1 xlO"* 


2SEc 


O AGB 


587 


O-rich 


RSG 


8.5 xlO* 


2SEc 


O AGB 


943 


O-rich 


RSG 


8.7 xlO* 


2SEc 


RSG 


1429 


O-rich 


RSG 


9.1 xlO" 


2SEc 


O AGB 


141 


O-rich 


RSG 


9.2 xlO"* 


2SEc 


RSG 


1330 


O-rich 


RSG 


9.3 xlO"* 


2SEc 


RSG 


897* 


O-rich 


RSG 


9.4 xlO* 


2SEc 


O AGB 


810"= 


O-rich 


RSG 


9.4 xlO*: 


2SEc 


RSG 


1204 


O-rich 


RSG 


1.0 xlO^ 


2SEc 


RSG 


264 


O-rich 


RSG 


1.1 xlO-"' 


2SEc 


RSG 


839 


O-rich 


RSG 


1.2 xlO^ 


2SEc 


RSG 


593 


O-rich 


RSG 


1.2 xlO^ 


2SEc 


RSG 


870 


O-rich 


RSG 


1.3 xlO^ 


2SEc 


RSG 


589* 


O-rich 


RSG 


1.3 xlO^ 


2SEc 


RSG 


1328 


O-rich 


RSG 


1.3 xlO^ 


2SEc 


RSG 


1189 


O-rich 


RSG 


1.3 xlO^ 


2SEc 


O AGB 


597 


O-rich 


RSG 


1.6 xlO^ 


2SEc 


RSG 


1150<* 


O-rich 


MW O AGB 


4.0 xlO'^ 


2SEa 


O AGB 


1686<* 


O-rich 


MW O AGB 


4.6 xlO'^ 


2SEa 


C AGB 


412^* 


O-rich 


M W O AGB 


6.0 xlO'^ 


2SEa 


O AGB 


1677'^ 


O-rich 


MW O AGB 


9.7 xlO^ 


2SEc 


O AGB 


807* 


O-rich 


OH/IR 


4.0 xlO"' 


3SB 


OH/IR 


1182 


O-rich 


OH/IR SG 


4.7 xIQS 


5? 




1326'= 


O-rich 


Peculiar 


5.5 xlO* 


5? 


PN 


890*= 


O-rich 


Peculiar 


7.1 xlO"* 


5? 


PN 


1384 


C-rich 


C AGB 


5.8 xlO^ 


3CE 


OH/IR 


1400 


C-rich 


C AGB 


7.1 xlO^ 


3C? 


OH/IR 


1007S 


C-rich 


C AGB 


8.2 xlO^ 


3CE 


OH/IR 


04374-6831'^ 


C-rich 


C AGB 


1.0 xlO* 


3CE 




1488 


C-rich 


C AGB 


1.1 xlO" 


3CE 


OH/IR 


438 


C-rich 


C AGB 


1.1 xlO" 


3CE 


OH/IR 


87 


C-rich 


C AGB 


1.1 xlO"* 


3CE 


OH/IR 


1120 


C-rich 


C AGB 


1.1 xlO^ 


3CE 


OH/IR 


1652 


C-rich 


C AGB 


1.2 xlO* 


3CR 


OH/IR 


218 


C-rich 


C AGB 


1.3 xlO* 


3CE 


OH/IR 
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Table 2 — Continued 



MSX 


Envelop G 


IRS 


^IR 


KSPW 


EVPOl 


NcimG 


rli pm 1 r V 


type 


UQ 




type 




[21 


f3l 
L^J 


[41 


[51 


[61 
L^J 


95 


C-rich 


C AGB 


1.4 xlO* 


3CE 


C IR 


1282 


C-rich 


C AGB 


1.5 xlO* 


3CE 


OH/IR 


220 


C-rich 


C AGB 


1.6 xlO* 


3CE 


C IR 


1592 


C-rich 


C AGB 


1.6 xlO* 


SCR 


_ 


1492 


C-rich 


C AGB 


1.9 xlO* 


3CE 


C IR 


775 


C-rich 


C AGB 


1.9 xlO"* 


C? 


C IR 


1794 


PAH 


Hii 


1.7 xlO"* 


5UE 


PN 


222 


PAH 


Hll 


3.8 xlO* 


5UE 


PN 


1306 


PAH 


Hii 


4.3 xlO* 


5UE 


PN 


894 


PAH 


Hii 


5.6 xlO* 


SUE 


PN 


22 


PAH 


Hii 


6.8 xlO* 


5UE 


PN 


836 


PAH 


Hii 


7.7 xlC* 


SUE 


PN 


217 


PAH 


Hii 


8.1 xlO* 


SUE 


PN 


559 


PAH 


Hii 


2.0 xlO^ 


SUE 


PN 


934 


PAH 


Hii 


2.5 xlO^ 


SUE 


PN 


764 


PAH 


Hii 


2.7 xlO^ 


SUE 


PN 


889 


PAH 


Hii 


4.0 Xl05 


SUE 


PN 


769' 










PN 


773' 










O AGB 



''For these objects, local sky subtraction was used. 

''Only SL module spectra were obtained for this object. The 
luminosity was calculated using MSX LMC 1189 as a template for 
the spectral shape, so is highly uncertain. 

'^This object exhibits extended emission due to a nearby or com- 
panion object (resolved in SL but not LL spectral images). Slit 
losses result in a large flux jump between the SL and LL modules, 
so the SL spectra were scaled by a factor of 1.86 to match the LL2 
spectrum; the luminosity is therefore very uncertain. 

'^The luminosities for these Milky Way O AGB stars were cal- 
culated using the following assumed distances: MSX LMC 1150 
4.2 kpc, MSX LMC 1686 3.3 kpc, MSX LMC 412 4.6 kpc, and 
MSX LMC 1677 5.1 kpc. These distances were estimated using 
the K-band distance modulus with absolute K-band magnitudes 
derived using the color-period-luminosity relations for Mira vari- 
ables (Whitelock et al. 1994). 

''The luminosities for these two objects were taJjen from Kastner 
et al. (2006). Ljpj^ obtained from the flux density integrated over 
the wavelength range 2.2 - 3S fJLVCi only. 

^A second source fell partially in the slit for several modules and 
nod positions of this observation. The extraction aperture used for 
this object was defined manually, using a tapering column of the 
default SMART widths to ensure accurate flux calibration. 
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^For the SLl module of this spectrum, only the second slit posi- 
tion was used in the final spectrum, as the flux density of the SLl 
spectrum in the first nod position did not match the rest of the 
modules and nod positions. 

''This object was not detected by MSX, so is designated by its 
IRAS name. 

'These two objects were not observed due to peak-up failure. 



Note. Columns: [1] MSX/IRAS name; [2] the dominant 
chemistry of the dust envelope, determined by the IR spectral fea- 
tures; [3] our classification of the likely type of object, based on 
the IRS spectral features, luminosity and all other available infor- 
mation; [4] the IR (1 - 100 /im) luminosity derived as described in 
§3.2. Uncertain luminosities are indicated by a colon; [5] the spec- 
tral class of the object in the scheme of Kraemer et al. (2002); [6] 
the type of the object determined from the NIR colors by EVPOl. 



-30- 



Table 3. Spitzer/ 2M ASS color diagnostics for object types. 



Class Color-color diagnostics 

H-[5.8] [8.0]-[24] K-[5.8] H-K 

[1] [2] [3] [4] [5] 



RSG 


-2.5 - 


-1.0 


AGB^ 


-0.5 - 


- 1.2 


MW AGB 


-2.3 - 


- 0.0 


C AGB 


0.5 - 


- 6.0 


HII 


2.0 - 


-4.5 



[ ^-[5.8] 
I 3 



+3.0] 



fH-\5.8] 



+0.5] 



3 +4.5] 
3.5 - 4.0 

-[5.8] 



L 3 
rg-[5.! 



+2.7] 
+2.5] 
4.5 - 8.0 



r if-JC-l.lS l 
L 0.3 J 



L 0.2 J 
r H-K-2.0 ] 
L 0.3 i 



H-K- 


-0.8j 


0.3 




-1.6 - 


0.0 


H-K- 


-0.9j 


0.2 




H-K- 


-1.2] 



0.3 J 

0.2 -3.0 



0.2 - 0.6 
0.7 - 1.3 
0.4 - 0.75 
1.1 - 2.7 
0.7 - 1.5 



Note. — Columns: [1] Type of IR object; [2]-[3] Color diagnostics for classifying objects based 
on Spitzer and 2MASS colors. 

^These criteria are based on the colors of only four objects and thus are very uncertain. 
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Fig. 1. — (a) Spitzer IRS spectra of the stars showing sihcate features. The spectra are 
arranged according to IRS spectral type and then object bolometric luminosity (Table 2). 
For each object, the left panel shows the IRS spectrum on a linear wavelength scale, and 
the right panel shows the spectrum with photometry from 2MASS, MSX, and IRAS, on a 
logarithmic wavelength scale. Alternating symbols are used to make clear which photometric 
points are associated with which IRS spectrum. The flux density is plotted on a logarithmic 
scale and the value in parentheses next to the object name indicates the offset in flux density 
applied to the spectrum. 



-32- 



7 - 



-1 



10 



20 30 
Wovelength (/xm) 



890 (+7.6) 



807 (+6.4) 




593 (+1.3) 



839 (+1.0) 
264 ( + 0.4) 



40 



1 6 - 



1 189 (+3.2) - 6 



14 - 



12 - 



10 ■ 



8 - 4--*- -A 



♦ 

i * A 




807 ( + 15.0) 



412 ( + 10.3) 
1686 (+9.2) 




10 

Wavelength (/xm) 



100 



Fig. 1. — (a) Continued 
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Fig. 1. — (b) Spitzcr IRS spectra of the stars showing SiC features. In the right panel, 
2MASS, MSX, and IRAS data are also shown where available. 
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Fig. 1. — (c) Spitzer IRS spectra of red objects showing emission features. In the right 
panel, 2MASS, MSX, and IRAS data are also shown where available. 
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Fig. 2. — A typical spectrum from each of the groups of objects in the sample, show- 
ing the primary spectral features detected for each, (a) The spectrum of 0-rich object 
MSX LMC 264. (b) The spectrum of C-rich object MSX LMC 1400. (c) The spectrum of 
emission-line object MSX LMC 22. The "jump" at ~14 //m is discussed in §4.3. 
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Fig. 3. — Figures illustrating the method of calculating the luminosities for the 0-rich (top), 
C-rich (middle) and emission-line (bottom) objects. The flux density was integrated under 
the IRS spectrum between 5 and 36 /im (solid line); a warm blackbody (dotted line), scaled 
to match the 2MASS photometry (open squares), was used to estimate the 1-5 /im flux 
densities; and a cool, emissivity- weighted blackbody (dashed line) was used to estimate the 
36 - 100 iim flux densities. 
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Fig. 4. — Histograms showing the distributions of the IR luminosity for the whole sample 
(upper panel) and for the 0-rich, C-rich and PAH-rich stars separately (lower panels). The 
hatched regions of the histogram for 0-rich objects indicate stars classified as RSGs based 
on their IR luminosities (diagonal). Stars identified as Galactic objects (see Table 2) are not 
shown. 
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Fig. 5. — Optical (DSS) and infrared (MSX A-band) images of the emission-line objects in 
the sample. Each image is 5' x 5'. The extended nebulous emission, particularly in the 
MSX 8.3 /im images, argues that these objects are H II regions associated with massive star 
formation regions, rather than planetary nebulae. 



-39- 




Fig. 6. — Comparison of the IRS spectrum of the emission-hne source MSX LMC 22 (middle) 
with the ISO spectra of PN NGC7027 (top), which has spectral subclass 4PU in the KSPW 
system, and Hll region IRAS10589-6034 (bottom), which has spectral class SUE. For clarity, 
the flux densities of the comparison sources have been scaled by the values shown. 
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Fig. 7. — 2MASS / MSX color-color diagram with symbols indicating the IRS spectral type: 
RSGs (open, inverted triangles), O AGB stars (open triangles), Galactic Mira variables (filled 
triangles), C AGB stars (open circles), H II regions (open squares), B stars (open stars), and 
OH/IR stars (open diamonds). The boxes indicate the new NIR color criteria to classify 
IR-luminous LMC objects. 
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Fig. 8. — Color-magnitude diagrams for our sample derived from 2MASS / MSX photometry. 
The symbols are the same as Figure 7 and indicate the IRS spectral type (listed in Table 
2). The filled symbols indicate Galactic objects. The dotted line in the right panel indicates 
the criterion to separate C AGB stars and H II regions. 
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Fig. 9. — Left: Spitzer synthetic IRAC / MIPS color-color diagram for the program objects, 
and some additional Galactic sources. The synthetic colors for the program objects were 
derived from the IRS spectra and the instrument spectral response functions. The symbols 
are the same as in Figure 7 and indicate the IRS spectral type (listed in Table 2). Additional 
(filled) symbols are colors derived from ISO spectra for Galactic 0-rich (filled triangles ) and 
C-rich (filled circles) stars, H II regions (filled squares), and PNe (filled stars) for compari- 
son. Also included is a line representing black bodies of different temperatures. Right: An 
expanded region of the left panel, showing a more limited color range and excluding the 
Galactic objects, to distinguish the regions occupied by the 0-rich RSGs and C-rich AGB 
stars more clearly. 
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Fig. 10. — Spitzer / 2MASS color-color diagrams showing the proposed criteria to classify 
IR-luminous LMC objects. Synthetic IRAC and MIPS fluxes were derived from the IRS 
spectra and the instrument spectral response functions. The symbols are the same as figure 
7 and indicate the IRS spectral type (listed in Table 2) . 
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Fig. 11. — Color- magnitude diagrams derived from synthetic Spitzer photometry. The sym- 
bols are the same as Figure 7. 
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Fig. 12. — The "bolometric correction" for the C-rich AGB stars in the sample, derived from 
the luminosity and K-band flux density, compared with the K-A color. The solid line shows 
a least-squares fit to the data points, and the dotted lines indicate the maximum error in 
the predicted luminosity (0.3 dex). 
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Fig. 13. — The "bolometric correction" for the RSGs in the sample, derived from the lumi- 
nosity and K-band flux density, compared with the K magnitude. The solid line shows the 
mean bolometric correction, and the dotted lines indicate the uncertainties determined from 
the sources with the largest distance from the mean. 



